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Abstract. In this paper, we present an optimal 0(n 2 ) time algorithm for deciding if a metric space (X,d) on n 
points can be isometrically embedded into the plane endowed with the ?i-metric. It improves the 0(n 2 log 2 n) 
time algorithm of J. Edmonds (2008). Together with some ingredients introduced by Edmonds, our algorithm 
uses the concept of tight span and the injectivity of the ii-plane. A different 0(n 2 ) time algorithm was recently 
proposed by D. Eppstein (2009). 

1. Introduction 

Deciding if a finite metric space (A", d) admits an isometric embedding or an embedding 
with a small distortion into a given geometric space (usually M fc endowed with some norm- 
metric) is a classical question in distance geometry which has some applications in theoretical 
computer science, visualization, and data analysis. The first question can be answered in 
polynomial time if M fc is endowed with the Euclidean metric due to classical results of Menger 
and Schonberg [6]. On the other hand, by a result of Frechet [6], any metric space can be 
isometrically embedded into some M fc with the ^-metric. However, it is NP-hard to decide 
if a metric space isometrically embeds into some M. k endowed with the l\ (alias rectilinear or 
Manhattan) metric [2,6]. More recently, Edmonds [9] established that it is even NP-hard to 
decide if a metric space embeds into M 3 with loo-metric (a similar question for M 3 with l\- 
metric is still open). In case of M 2 , both l\- and Zoo-metrics are equivalent because the second 
metric can be obtained from the first one by a rotation of the plane by 45° and then by a 
shrink by a factor -4=. The embedding problem for the rectilinear plane was investigated in 
the papers [3,12], which ultimately show that a metric space (X, d) embeds into the Zi-plane 
if and only if any subspace with at most six points does [3] (a similar result for embedding 
into the Zi-grid was obtained in [4]). As a consequence, it is possible to decide in polynomial 
time if a finite metric space embeds into the Zi-plane. Edmonds [9] presented an 0(n 2 log 2 n) 
time algorithm for this problem and very recently we learned that Eppstein [10] described 
an optimal 0(n 2 ) time algorithm (for earlier algorithmic results, see also [5]). In this note, 
independently of [10], we describe a simple and optimal algorithm for this embedding problem, 
which is different from that of [10]. 

We conclude this introductory section with a few definitions. In the sequel, we will denote 
by di or by || • ||i the Zi-metric and by doc the £oo- m etric. A metric space (X, d) is isometri- 
cally embeddable into a host metric space (Y, d') if there exists a map <p : X — > Y such that 



d'(ip(x), (f(y)) = d(x, y) for all x, y € X. In this case we say that X is a subspace of Y. A retrac- 
tion ip of a metric space (Y, d) is an idempotent nonexpansive mapping of Y into itself, that 
is, ip 2 = ip : Y — > Y with d(tp(x),<p(y)) < d(x,y) for all x,y G y. The subspace of Y induced 
by the image of Y under ip is referred to as a retract of Y. Let (X, d) be a metric space. The 
(closed) ball and the sphere of center x and radius r are the sets B(x, r) = {p £ X : d(x,p) < r} 
and S(x,r) = {p £ I : d(x,p) = r}, respectively. The interval between two points x,y of X is 
the set I(x,y) = {z £ X : d(x,y) = d(x,z) +d(z,y)}. Any ball of (R fe , d^) is an axis-parallel 
cube. A subset S of X is gated if for every point x £ X there exists a (unique) point x' £ S 1 , 
the yaie of x in 5, such that x' € y) for all y € S [8]. The intersection of gated sets is also 
gated. A geodesic in a metric space is the isometric image of a line segment. A metric space is 
called geodesic (or Menger- convex) if any two points are the endpoints of a geodesic. 

For a point p of M 2 , denote by Qi(p), ■ ■ ■ , Qi(p) the four quadrants of M 2 defined by the 
vertical and horizontal lines passing via the point p and labeled counterclockwise. Any interval 
I\(x,y) of the rectilinear plane (M 2 ,di) is an axis-parallel rectangle which can be reduced to a 
horizontal or vertical segment. Any ball of (M 2 ,di) is a lozenge obtained from an axis-parallel 
square by a rotation by 45° degrees. In the rectilinear plane, any halfplane defined by a vertical 
or a horizontal line is gated. As a consequence, axis-parallel rectangles, quadrants, and strips 
of (M 2 ,c?i) are gated as intersections of such halfplanes. 

2. Tight spans 

A metric space (X,d) is called hyperconvex (or injective) [1,11] if any family of closed balls 
B(xi,ri) with centers Xj and radii rj, i 6 I, satisfying d(xi,Xj) < rj + rj for all i,j G I has a 
nonempty intersection, that is, (X, d) is a geodesic space such that the closed balls have the 
Helly property. Since the closed balls of (M fc ,doo) are axis-parallel boxes, the metric spaces 
(R fc ,doo) an d (M. 2 ,d\) are hyperconvex. It is well known [1] that (X,d) is hyperconvex iff it 
is an absolute retract, that is, (X, d) is a retract of every metric space into which it embeds 
isometrically. As shown by Isbell [11] and Dress [7], for every metric space (X,d) there exists 
the smallest injective space T(X) extending (X,d), referred to as the injective hull [11], or tight 
span [7] of (X,d). The tight span of a finite metric space (X,d) can be defined as follows. Let 
T(X) be the set of functions / from X to M such that 

(1) for any x,y in X, f(x) + f(y) > d(x,y), and 

(2) for each x in X, there exists y in X such that f(x) + /(y) = d(x,y). 

One can interpret f{x) as the distance from / to x. Then (1) is just the triangle inequality. 
Taking x = y in (1), we infer that fix) > for all x G X. The requirement (2) states 
that T(X) is minimal, in the sense that no value f(x) can be reduced without violating the 
triangle inequality. We can endow T(X) with the Zoo-distance: given two functions / and g 
in T(X), define p(f,g) = max|/(x) — g(x)\. The resulting metric space (T(X),p) is injective 
and (T(X),p) is called the tight span of (X,d). There is an isometric embedding of X into its 
tight span T(X). Moreover, any isometric embedding of (X,d) into an injective metric space 
(Y,d') can be extended to an isometric embedding of (T(X), p) into (Y,d'), i.e., (T(X),p) is the 
smallest injective space into which (X, d) embeds isometrically. 

In general, tight spans are hard to visualize. Nevertheless, if \X\ < 5, Dress [7] completely 
described T(X) via the interpoint-distances of X. For example, if \X\ = 3, say X = {x,y, z}, 
then T(X) consists of three line segments joined at a (Steiner) point, with the points of X 
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Figure 2. The three canonical types of tight span of 5-point metric space. 



at the ends of the arms (see Fig. CD (a)). The lengths of these segments are a x ,a y ,a z , where 
Oix '■= {y, z) x = l/2(d(x, y) + d(x, z) — d(y, z)) is the Gromov product of x with the couple y, z 
[a y and a z are defined in a similar way). Notice that one of the values a x ,a y , a z may be 0, in 
this case one point is located between two others. If \X\ = 4, then the generic form of T{X) 
is a rectangle R{X) endowed with the Zi-metric, together with a line segment attached by one 
end to each corner of this rectangle (see Fig. [T] (b)). The four points of X are the outer ends of 
these segments. The lengths of these segments and the sides of the rectangle can be computed 
in constant time from the pairwise distances between the points of X; for exact calculations 
see [7]. It may happen that R(X) degenerates into a segment or a point. Finally, there are 
three canonical types of tight spans of 5-point metric spaces precisely described in [7] (see also 
Fig. [2] for an illustration). Each of them consists of four or five rectangles, five segments, and 
eventually one rectangular triangle, alltogether constituting a 2-dimensional cell complex. All 
sides of the cells can be computed in constant time as described in [7]. It was also noticed 
in [7] that if for each quadruplet X' of a finite metric space (X, d) the rectangle R(X') is 
degenerated, then (X, d) isometrically embeds into a (weighted) tree and its tight span T(X) 
is a tree-network. 

From the construction of tight spans of 3- and 4-point metric spaces immediately follows that 
any metric space (X, d) with at most 4 points and its tight span (T(X), p) can be isometrically 
embedded into the Zi-plane as shown in Fig. [T](b). This is no longer true for metric spaces on 5 
points: to embed, some cells of the tight span must be degenerated. If \X\ = 4 and the rectangle 
R(X) is non-degenerated, one can easily show that R{X) isometrically embeds into the Zi-plane 
only as an axis-parallel rectangle. Therefore, if additionally the four line segments of T{X) are 
also non-degenerated, then up to a rotation of the plane by 90°, X and T(X) admit exactly 
two isometric embeddings into the /i-plane. If one corner of R(X) is a point of X and the 
embedding of the rectangle R(X) is fixed, then there exist three types of isometric embeddings 
of X and T(X) into the rectilinear plane: two segments of T(X) are embedded as axis-parallel 
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segments and the third one as a segment whose slope has to be determined. Analogously, if two 
incident corners of R(X) are points of A, the two segments of T(X) are either embedded as 
axis-parallel segments, or one as a horizontal or vertical segment and another one as segment 
whose slope has to be determined. Note also that from the combinatorial characterization of 
finite metric subspaces of the Zi-plane presented in [3] immediately follows that a tree-metric 
(X, d) is isometrically embeddable into the /i-plane if and only if the tree-network T(X) has at 
most four leaves. Finally note that since (K 2 ,<ii) is injective, by minimality property of tight 
spans, T(X) is an isometric subspace of the /i-plane for any finite subspace X of M 2 . 

3. Algorithm and its correctness 

3.1. Outline of the algorithm. Let (X, d) be a metric space with n points, called terminals. 
Set X = {x\, . . . ,x n }. Our algorithm first finds in 0(n 2 ) time a quadruplet P° of X whose 
tight span contains a nondegenerated rectangle R(P°). If such a quadruplet does not exists, 
then (X, d) is a tree-metric and T{X) is a tree-network. If this tree-network contains more 
than four leaves, then (A, d) cannot be isometrically embedded into the /i-plane, otherwise 
such an embedding can be easily derived. Given the required quadruplet P°, we consider any 
isometric embedding of P° and of its tight span into the Zi-plane as illustrated in Fig. |4] and 
partition the remaining points of A into groups depending on their location in the regions of 
the plane defined by the rectangle R(P°) and the segments of T(P°). The exact location of 
points of A in these regions is uniquely determined except the four quadrants defined by R(P°). 
At the second stage, we replace the quadruplet P° by another quadruplet P by picking one 
furthest from R(P°) point of A in each of these quadrants. We show that the rectangle R(P°) 
is contained in the rectangle R(P), moreover, for any isometric embedding ipo of P and T{P) 
into the li-plane, the quadrants defined by two opposite corners are empty (do not contains 
other terminals of A). Again the location of the points of A in all regions of the plane except 
the two opposite quadrants is uniquely determined. To compute the location of the remaining 
terminals in these two quadrants we adapt the second part of the algorithm of Edmonds [9]: 
we construct on these terminals a graph as in [9], partition it into connected components, 
separately determine the location of the points of each component, and then combine them 
into a single chain of components in order to obtain a global isometric embedding ip of (A, d) 
extending (£>o or to decide that it does not exist. 

Now, we briefly overview the algorithms of Edmonds [9] and Eppstein [10]. Edmonds [9] 
starts by picking two diametral points p, q of A. These two points can be embedded into the l±- 
plane in an infinite number of different ways. Each embedding defines an axis-parallel rectangle 
IT whose half-perimeter is exactly d(p, q). Using the distances of p and q to the remaining points 
of A, Edmonds computes a list A of linear size of possible values of the sides of the rectangle 
II. For each value 5 from this list, the algorithm of [9] decides in 0(n 2 ) time if there exists 
an isometric embedding of A such that one side of the rectangle II has length 5. For this, it 
partitions the points of A into groups, depending on their location in the regions of the plane 
determined by IT. In order to fix the positions of points in one of these regions, Edmonds [9] 
defines a graph whose connected components are also used in our algorithm. While sweeping 
through the list A, the algorithm of [9] update this graph and its connected components in 
an efficient way. Notice that the second part of our algorithm is similar to that from [9], 
but instead of trying several sizes of the rectangle II, we use the tight spans to provide us 
with a single rectangle, ensuring some rigidity in the embedding of the remaining points. The 
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Figure 3. The tree- network Aj,. 

algorithm of Eppstein [10] is quite different in spirit from our algorithm and that of Edmonds [9]. 
Eppstein [10] first incrementally constructs in 0(n 2 ) time a planar rectangular complex which 
is the tight span of the input metric space (X, d) or decide that the tight span of X is not 
planar. In the second stage of the algorithm, he decides in 0(n 2 ) time if this planar rectangular 
complex can be isometrically embedded into the Zi-plane or not. 

3.2. Computing the quadruplet P° . For each i = 1, . . . , n, set Xi := {xi, . . . , x{\. We start 
by computing the tight span of the first four points of X. If this tight span is not degenerated 
then we return the quadruplet X\ as P°. Now suppose that the tight span of the first i — 1 
points of X is a tree-network A^i with at most four leaves. This means that ^4j_i contains 
one or two ramification points (which are not necessarily points of X) having degree at most 
4, all remaining terminals of Xi—i are either leaves or vertices of degree two of Ai—\. We say 
that two terminals of are consecutive in A\-\ if the segment connecting them in ^4i_i 

does not contain other points of Note that ^4i_i contains at most n + 4 of consecutive 

pairs. For each pair Xj,x k of consecutive terminals of we compute the Gromov product 
a Xi := (xj,Xk) Xi = l/2(d(xi,Xj) + d(xi,Xk) -d(xj,x k )) of Xi with {xj,x k }. Let {a,b} be the 
pair of consecutive points of minimizing the Gromov product a Xi = (a, b) Xi . Let c be the 
point of the segment [a, b] of Ai_i located at distance a a := (b,Xi) a from a and at distance 
ctb := (a,Xi)b from b (c may coincide with one of the points a or b). 

Denote by Ai the tree-network obtained from Ai_\ by adding the segment [xi,c\ of length 
a Xi . By running Breadth-First-Search on Ai rooted at x^, we check if dA t { x ii x j) = d(xi,Xj) for 
any terminal Xj of Xi . If this holds for all Xj € Xi , then the tight span of Xi is the tree- network 
Ai. If Ai contains more than 4 leaves, then we return the answer "not" and the algorithm halts. 
Otherwise, if i = n, then we return the answer "yes" and an isometric embedding of X and 
its tight span A n in the ^i-plane, else, if i < n, we consider the next point Xi + \. Finally, if 
Xj is the first point of Xi such that dAi(xi,Xj) ^ d(xi,Xj), then we assert that the tight span 
of the quadruplet {a,b,Xi,Xj} is non- degenerated and we return it as P°. Suppose by way of 
contradiction that T(P°) is a tree. Since Ai_i realizes and T(P°) realizes P°, the subtree 
of A^i spanned by the terminals a, b, and Xj is isometric to the subtree of T(P°) spanned by 
the same terminals. On the other hand, T(P°^ contains a point c' located at distance Q^,cy a , 
and at, from Xi,a, and b, respectively. This means that T(P°) is isometric to the subtree 
of Ai spanned by the vertices Xi,a,b, and Xj, (see Fig. [3|) contrary to the assumption that 
dAi(xi,Xj) ^ d(xi,Xj). Hence, this inequality implies indeed that T(P°) is not a tree. Finally 
note that dealing with a current point Xi takes time linear in i, thus the whole algorithm for 
computing the quadruplet P° runs in 0(n 2 ) time. 
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3.3. Classification of the points of X with respect to the rectangle of T(P°). Let 
P° = {pl,p2,pl,p^} be the quadruplet whose tight span T{P°) is non-degenerated. Let R° be 
one of the two possible isometric embeddings of the rectangle R(P°) of T(P°) and consider 
a complete or a partial isometric embedding of T(P°) such that R(P°) is embedded as R°. 
Denote by Q\ , Q\ , Q\ , Q\ the four (closed) quadrants defined by the four consecutive corners 
?i, 92> ?3> ?4 °f ^° labeled in such a way that the point p° must be located in the quadrant 
Q°,i = 1, ...,4. Let also , S 1 ^ , £3 , and 5| be the remaining half-infinite strips. Since we 
know how to construct in constant time the tight span of a 5-point metric space, we can 
compute the distances from all terminals p of X to the corners of the rectangle R(P°) (and 
hence to the corners of R°) in total 0{n) time. With some abuse of notation, we will denote 
the ^-distance from p to the corner q° of R° by d(p, q°). Since R° is gated, from the distances 
of p to the corners of R° we can compute the gate of p in R°. Consequently, for each point 
p £ X \ P° we can decide in which of the nine regions of the plane will belong its location (p(p) 
under any isometric embedding (p of (X, d) subject to the assumption that R(P°) is embedded 
as R°. If <p(p) belongs to one of the four half-strips or to R°, then we can also easily find 
the exact location itself: this can be done by using either the gate of p in R° or the fact that 
inside these five regions the intersection of the four /i-spheres centered at the corners of R° 
and having the distances from respective corners to p as radii is a single point. So, it remains 
to decide the locations of points assigned to the four quadrants Q\-,Q 21 Q% an d Q\- For any 
point p E X which must be located in the quadrant Q°, the set of possible locations of p is 
either empty (and no isometric embedding exists) or a segment s p of Q° consisting of all points 
z e Q° such that \\z - q°\\i = d(p, q°). 

Notice that for any quadruplet P' = {p'i^P^^v'^^p'a} of terminals such that p[ is assigned to 
the quadrant Q°, i € {1,2,3,4}, the rectangle R° belongs to the tight span T(P') of P' . Indeed, 
for any point p[, i G {1, 2, 3, 4} and any point r of R°, we have — r||i + ||r — p'j\\i = — p'j\\i, 
where j is selected in such a way that q° and q? are opposite corners of R°. From injectivity of 
the ?i-plane and the characterization of tight spans we conclude that all points of R° belong 
to T(P'), establishing in particular that this tight span is also non-degenerated. 

3.4. The quadruplet P and its properties. Let P = {pi,P2)P3>f?4} be the quadruplet of 
X, where pi is a point of X which must be located in the quadrant Q° and is maximally distant 
from the corner q° of R°. As we established above, the tight span of P is non-degenerated, 
moreover the rectangle R(P) contains the rectangle R(P°). As we also noticed, there exists a 
constant number of ways in which we can isometrically embed T(P) into the ^i-plane. Further 
we proceed in the following way: we pick an arbitrary isometric embedding cpo of T(P) and try 
to extend it to an isometric embedding <p of the whole metric space (X, d) in the /i-plane. If 
this is possible for some embedding of T(P), then the algorithm returns the answer "yes" and 
an isometric embedding of X, otherwise the algorithm returns the answer "not" . Let R denote 
the image of R(P) under ifQ. 

We call a terminal pi of P fixed by the embedding (fa if either fo(pi) is a corner of the 
rectangle R or the segment of T(P) incident to pi is embedded by (po as a horizontal or a 
vertical segment; else we call pi free. The embedding of a free terminal pi is not exactly 
determined but is restricted to a segment s Pi consisting of the points of the quadrant defined 
by qi and having the same Zi-distance to We call the terminals Pi,Pj+i( mo d4) incident and 
the terminals Pi,Pi+2(mod4) opposite. From the isometric embedding of T(P) we conclude that 
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Figure 4. Possible isometric embeddings of T(P). 
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Figure 5. The partition of the plane into half-strips and quadrants. 



at most one of two incident terminals can be free. Moreover, if a terminal pi of P is fixed but is 
not a corner of R, then at least one of the two terminals incident to pi is also fixed. If all four 
tips of T(P) are non-degenerated, then all four terminals of P are fixed. If only three tips of 
T(P) are non-degenerated then at most one terminal of P is free, all remaining terminals are 
fixed. If only two tips of T{P) are non-degenerated, then either they correspond to incident 
terminals, one of which is fixed and another one is free or to two opposite terminals which are 
both free. Finally, if only one tip of the tight span is non-degenerated, then it corresponds to 
a free terminal, all other terminals of P are corners of R and therefore are fixed (see Fig. 0] for 
the occurring possibilities). 

Denote by II the smallest axis-parallel rectangle containing R and the fixed terminals of 
P; Fig. [5] illustrates II for two cases from Fig. H] (if a terminal is free, then the respective 
corner of R is also a corner of II). Let qi,q2,Q3,Q4 be the corners of II labeled in such a way 
that qi is the corner of R corresponding to the point pi and to the corner q° of R°. Denote by 
Qi, ■ ■ ■ ,Qi the quadrants of E? defined by the corners of II and by S\, . . . ,Si the remaining 
half-infinite strips. Again, as in the case of the quadruplet P°, by building the tight spans of 
P U {p} for all terminals p £ X \ P, we can compute in total linear time the distances from 
all such points p to the corners of R (and to the corners of II). From these four distances 
and the distances of p to the terminals of the quadruplet P we can determine in which of the 
nine regions Qi, Q2, Q3, Q4, Si, S2, S3, £4, II of the plane must be located p. Moreover, if p is 
assigned to the rectangle II or to one of the four half-strips Si , S2 , S3 , S4, then we can conclude 
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Figure 6. On possible locations of terminals in Q\ and Q3. 

that, in the region in which p assigned, the intersection of the four spheres centered at the 
terminals of P and having the distances from respective points to p as radii is either empty or 
a single point. The sphere centered at a free terminal pi is needed only to decide the location 
of p in the quadrant Q of the plane having the same apex a as the quadrant Q, t and which is 
opposite to Qi (a is a corner of II). But in this case, instead of considering the sphere of radius 
d{p, Pi) centered at po(pi) we consider the sphere of radius d(p,pi) — ||y»o (Pi) — and centered 
at a : indeed, both these spheres have the same intersection with Q. 

We are now ready to prove the following property of the quadruplet P : among the four 
quadrants Qi,Q2,Q3, o.nd Q4 defined by P, two opposite quadrants are empty, i.e., they do not 
contain terminals of X\P. First note that by inspecting the different cases listed in Fig. H] one 
can check that the two neighbors Pj_i( mo d4) an d Pi+i(mod4) °f a f ree point pi G P are both fixed; 
let say p\ and pj, are fixed. Now, suppose by way of contradiction that a terminal q £ X \ P 
must be located in the quadrant Q\. This means that its gate in the rectangle II is the corner 
of IT corresponding to p\. Since in any embedding ip of X that extends the chosen embedding 
of T(P°) the terminal p\ is located in Q°, we deduce that Qi(<p(pi)) Q Q°- On the other hand, 
the inclusion Q\ C Qi{p{p\)) follows directly from the definition of Q\ and the fact that p\ is 
fixed. Now, from the inclusions Q\ C Qi(tp(pi)) C QJ, we obtain that q € Q\ and, since q is 
closer to p\ than to q\ , we get a contradiction with the choice of p\ , establishing that indeed 
Qi does not contain any point of X \ P. The same argument shows that Q3 is empty as well. 
Note that actually we proved that any quadrant Qi corresponding to a fixed terminal pi of P 
is empty. 

3.5. Locating in the non-empty quadrants Q\ and Q3. As we have showed in previous 
subsection, any isometric embedding (p of (X, d) extending the embedding tpo of T(P) locates 
each terminal p of X \ P in one and the same of the nine regions defined by II. Moreover, if p 
must be located in the rectangle II or in one of the four half-strips Si, ■ ■ ■ , 64, then this location 
<p(j)) is uniquely determined from the distances to the terminals of P and to the corners of II. 
We also established that only one or two opposite quadrants defined by II, say Q\ and Q3, can 
host terminals of X \ P; see Fig. [6l We will show now how to find the exact location of the set 
X\ of terminals assigned to Q\ (the set X3 of terminals which must be located in Q3 is treated 
analogously) . 

Note that independently of how the extension ip of (po is chosen, for each terminal u € X\, 
the /i-distance ||y?(it) — from the location of u to the corner q\ of II is one and the same, 
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Figure 7. tp(u) and ip(v) are 
fixed by tp(p*) and <p(p*) 



Figure 8. ip(v) is fixed by <p{p*) 



which we denote by A u . The value of A u can be easily computed because q\ lies between (p(u) 
and ip(pi) for any pi £ P: for example, we can set A u := d(u,pi) — ||</?o(Pi) ~~ <7i||i- Then the 
set of all possible locations (f(u) of it 6 Xi is the Zeue/ segment s u which is the intersection of 
Qi with the sphere S(q\, A u ) of radius A u centered at q\. 

To compute the locations of the terminals of X\ in the quadrant Qi, we adapt to the Zi-plane 
the definition of a graph (which we denote by G\ = (X\,Ei)) defined by Edmonds [9] in the 
Zoo-plane. Two terminals it, i; £ X\ are adjacent in G\ if and only if d(u,v) > \A U — A v \. 
Equivalently it, v £ X\ with A u < A v are adjacent in G\ iff it cannot be located between q\ 
and v : ip(u) £ Ii(qi, <f(v)). Denote by C±, C2, . . . ,Ck the connected components of the graph 
G\. They have the following useful properties established in Lemmata 3-5 of [9]: 

(1) Each component Cj is rigid, i.e., once the location of any point u of Q has been fixed, the 
locations of the remaining points of Cj are also fixed (up to symmetry with respect to the line 
parallel to the bisector of Q\ and passing via it); 

(2) The components C\, C2, ■ ■ ■ , of the graph G\ can be numbered so that the points of each 
Ci appear consecutively in the list of points it € X\ sorted in increasing order of their distances 
A u to 91 ; 

(3) For a component Cj of G±, let Bi be the smallest axis-parallel rectangle containing {(fi(u) : 
it G Ci} for an isometric embedding ipi of (Ci,d) in the Zi-plane. Let b{ be the upper right 
corner of B{. Then the embedding of C\, C2, ■ ■ ■ , C/. preserves the distances between all pairs 
of points lying in different components if and only if for every pair of consecutive components 
Ci and Cj+i, the rectangle -Bj+i lies entirely in the quadrant Qiipi). 

The location in the quadrant Q\ of some terminals of X\ (and therefore of the connected 
components containing them) can be fixed by terminals already located in the two half-strips 
incident to Q\. We say that a terminal it € X\ is fixed by a terminal p already located in £iUS4 
if the intersection of the segment s u with the sphere S((p(p), d(p, it)) is a single point. Note that 
if it £ X\ is fixed by a terminal located in S\, then it is also fixed by the upmost terminal p* 
located in this half-strip. Analogously, if it £ X\ is fixed by a terminal of £4, then it is also fixed 
by the rightmost terminal p* located in S4. Therefore by considering the intersections of the 
segments s u ,u £ Xi, with the spheres S(<p(p*),d(p*,u)) and S (<p(p*) , d(p* , u)) we can decide 
in linear time which terminals of X\ are fixed by p* and and find their location in Q\ (for 
an illustration, see Fig. [7]). According to property (1), if a terminal of a connected component 




Figure 9. On the assemblage of blocks -Bj+i, . . . , B^. 



of G\ is fixed, then the location of the whole component is also fixed (up to symmetry). Let 
Cj be the connected component of G\ containing the furthest from q\ terminal u 6 X\ fixed 
by p* or p*, say by p* (therefore the location of Cj is fixed). We assert that all terminals of 
C\, C2, Cj-i are also fixed by p* . Indeed, pick such a terminal v. From property (2) we 
conclude that A v < A u and from the definition of G\ we deduce that v must be located in 
the axis-parallel rectangle I\{q\, <p(u)), and therefore below u. Since u is fixed by p* , u must 
be located below p*, whence v also must be located below p*. We can easily see that the 
intersection of s v with the sphere S((p(p*), d(p* ,v)) is a single point, i.e. v is also fixed by p* 
(see Fig. E}. 

It remains to locate in Q\ the terminals of the components Cj+i, Cj+2j • • • , We compute 
separately an isometric embedding of each component Ci for i = j + 1, . . . , k. For this, we fix 
arbitrarily the location of the first two points u, v of Ci in the segments s u and s v so that to 
preserve the distance d(u,v) (the terminals of Ci are ordered by their distances to q±). By 
property (1) of [9], the location of the remaining points of Ci is uniquely determined and each 
point w of Ci will be located in its level segment s w . Let ipi be the resulting embedding of 
Cj. Denote by Bi the smallest axis-parallel rectangle (alias box) containing the image <fi(Ci) 
of Ci. Let a, and bi denote the lower left and the upper right corners of Bi. Note that 
belongs to the ^-interval between q\ and the image y>i{u) of any terminal u of Cj, while the 
Zi-interval between gi and bi will contain the images of all terminals of Cj. Therefore if we set 
A ai := A u — \\a,i — (fi(u) ||i and := A u + \\tp(u) — bi ||i, where u is any terminal of Q, then in 
all isometric embeddings of (Ci,d) in which all terminals u € Ci are located on s u , the points 
aj and bi must be located on the level segments s ai and s^ , defined as the intersections of the 
quadrant Q\ with the spheres S(qi,A ai ) and S(qi, A^). 

By properties (2) and (3) of [9], in order to define a single isometric embedding of the 
components Cj+i, . . . , we now need to assemble the boxes -Bj+i, . . . , B^ (by moving their 
terminals along the level segments) in such a way that for two consecutive components Ci and 
Cj+i, the box Bi + \ lies entirely in the quadrant Q\(bi). We assert that this is possible if and only 
if for each pair of consecutive boxes Bi,Bi + i, i = j,j + 1, k — 1, the inequality A&. < A ai+1 
holds. Indeed, if A^ < A ai+1 , then translating Bi + \ along the segment s ai+1 , we can locate its 
corner cii+i in the quadrant Q\{bi) and thus satisfy the embedding requirement. Conversely, 
if A;,. > A ai+1 holds, then aj+i cannot belong to the quadrant Q\(pi) independently of the 
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positions of a.; + i and bi on their level segments. This local condition depends only of the values 
of A Qi , and is independent of the actual location of the boxes Bi, i = 1, . . . , k. As a result, 
the algorithm that embeds the boxes Bj + i, ...,-Bfc is very simple. For each i = j, k — 1, we 
compute the box -B,+i and the values of A ai+1 and A{>- +1 . If A ai+1 < A^ for some i, then return 
the answer "there is no isometric embedding of {X, d) extending the embedding ipo of T(P)". 
Otherwise, having already located the box Bi, by what has been shown above, the intersection 
of the quadrant Q\{bi) with the level segment s ai+1 is non-empty. Therefore we can translate 
Bj + \ in such a way that its lower left corner aj + i becomes a point of this intersection. 

In this way, we obtain an embedding of Cj+i , . . . , Ck and -Bj+i, . . . , Bf~ satisfying the con- 
ditions (l)-(3), thus an isometric embedding of the metric space (U(Lj+i Ci, d) in Q\. Anal- 
ogously, by constructing the graph G3 = (A3, £3) and its components, either we obtain a 
negative answer or we return an isometric embedding of the metric space defined by the non- 
fixed components of G3 in the quadrant Q3. Denote by (p the embedding of X which coincides 
with ifo on P, with these two embeddings on the non-fixed components of Gi and G3, and 
with the already computed fixed locations of the terminals assigned to II, to the half-strips 
Si, 1S2, S3, S4, and to the fixed connected components of the graphs Gi and G3. In 0(n 2 ) we 
test if (p is an isometric embedding of (X, d) into the Zi-plane. If the answer is negative, then 
we return "there is no isometric embedding of (X, d) extending the embedding ip$ of T(P)", 
otherwise we return </? as an isometric embedding. The algorithm returns the global answer 
"not" if for all possible embeddings ipo of T(P) it returns the negative answer. From what we 
established follows that in this case {X, d) is not isometrically embeddable into the /1-plane. 

3.6. Algorithm and its complexity. We conclude the paper with a description of the main 
steps of the algorithm and their complexity. 

Algorithm Embedding into the Zi-plane 

Input: A metric space (X, d) on n points 

Output: An isometric embedding if of (X, d) into (R 2 ,di) or the answer "not" if it does not exist 
Step 1. Find a quadruplet P° of X whose tight span contains a rectangle. If P° does not exist, then T(X) is 
a tree. If T(X) has more than 4 leaves, then return "not", else return an embedding of T(X) and (X, d). 
Step 2. Pick any embedding of T(P°) and for each terminal of X \ P° determine in which of the nine regions 
of the plane it must be located. Using this partition of X \ P° , define the quadruplet P. 

Step 3. Embed P and its tight span T(P) into the Zi-plane in all possible different ways. Try to extend each of 
these embeddings to an isometric embedding of (X, d) following the rules (a)-(g). If all of these attempts return 
the answer "not" , then return the answer "not" , else return one of the obtained embeddings. 

(a) Given an embedding ipo of T(P), for each terminal u of X \ P determine in which of the nine regions 
defined by the rectangle II will be located u in any isometric embedding extending ipo', 

(b) Locate the terminals assigned to the rectangle LT and the four half-strips Si, S2, S3, S4; 

(c) Define the sets of terminals X\ and X3 assigned to the quadrants Q\ and Q3, construct the graphs 
Gi — (Xi,Ei) and G3 = {Xz,Ez) and their connected components; 

(d) Find the terminals of X\ fixed by p* ,p* and their location in Q\. Do a similar thing for X3; 

(e) Find an isometric embedding of each component d of Gi not containing already fixed terminals so 
that its terminals are located on their level segments. Do a similar thing for G3; 

(f) Test if the free components Gj+i , . . . , Ck of Gi satisfy the condition < A ai+1 for i = j + 1, . . . , k — 1. 
If not, then return the answer "not", else locate consecutively the boxes Bj+i, . . . , in such a way 
that eii+i is located in Qi(bi) ns a;+1 and fix in this way the position of all terminals of X\. Do a similar 
thing for the free components of G3; 

(g) Verify if the resulting embedding of X extending ipo is an isometric embedding of (X, d). If "yes", 
then return it as a resulting isometric embedding, otherwise return the answer "there is no isometric 
embedding of {X,d) extending the embedding ipo" ■ 
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In Subsection 13,21 we established that the quadruplet P°, if it exists, can be computed in 
<9(n 2 ) time. If P° does not exists, then the tree-network A n (constructed within the same time 
bounds) is the tight span of (X,d). Embedding A n (if it has at most 4 leaves) in the Zi-plane 
can be easily done in linear time. As shown in Subsection 13.31 Step 2 can be implemented 
in linear time. There exists a constant number of ways in which the quadruplet P and its 
tight span can be isometrically embedded in the Zi-plane. Therefore, to show that Step 3 has 
complexity 0(n 2 ), it suffices to estimate the total complexity of the steps (a)-(g) for a fixed 
embedding <po of T{P). Step (a) is similar to Step 2, thus its complexity is linear. The exact 
location of each terminal in the half-strips or in II is determined as the intersection of two 
spheres, therefore step (b) is also linear. Defining the graph G\ and computing its connected 
components can be done in 0(|Xi| 2 ) time. Thus step (c) has complexity <9(n 2 ). Steps (d) and 
(e) can be implemented in an analogous way as (b), thus their complexity is 0(n). Testing 
the condition in step (f) and assembling the free components into a single chain is linear as 
well. Finally, step (g) requires 0(n 2 ) time. Therefore, the total complexity of the algorithm is 
<3(n 2 ). Summarizing, here is the main result of this note: 

Theorem 1. For a metric space (X,d) on n points, it is possible to decide in optimal 0(n 2 ) 
time if (X,d) is isometrically embeddable into the l\-plane and to find such an embedding if it 
exists. 
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